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Abstract. This work describes a method for esti- 
mating vertical fluxes of horizontal momentum carried 
by short horizontal scale gravity waves (X x = 10-100 
km) using aircraft measured winds in the lower strato- 
sphere. We utilize the in situ wind vector and pressure 
altitude measurements provided by the Meteorological 
Measurement System (MMS) on board the ER-2 air- 
craft to compute the momentum flux vectors at the 
flight level above deep convection during the tropical 
experiment of the Stratosphere Troposphere Exchange 
Project (STEP-Tropical). Data from Flight 9 are pre- 
sented here for illustration. The vertical flux of horizon- 
tal momentum in these observations points in opposite 
directions on either side of the location of a strong con- 
vective updraft in the cloud shield. This property of in- 
ternal gravity waves propagating from a central source 
compares favorably with previously described model re- 
sults. 

1. Introduction 

The evolving picture of how tropospheric wave sources 
control the structure and dynamics of the middle atmo- 
sphere now includes the important role played by high 
frequency gravity waves. Although typical power spec- 
tra of horizontal velocity or temperature perturbations 
are red in frequency, suggesting that the low frequency 
waves hold most of the energy, the high frequency waves 
carry a disproportionate share of the gravity wave mo- 
mentum flux. This is particularly true in the meso- 
sphere, where Fritts and Vincent [1987] found two- 
thirds of the momentum flux was carried by the high- 
est frequency band they observed: waves with periods 
between 8 minutes and an hour. It is the important 
quantity of momentum flux (and its divergence) that 
describes the mean flow acceleration due to these waves, 
which may be driving the middle atmosphere away from 
its radiatively determined state, and may play a role in 
forcing the semiannual oscillation (SAO) in the tropical 
mesosphere. (See e.g, Andrews et al. [1987] for details 
and references.) There is in general a lack of informa- 
tion about these shorter gravity waves which are rela- 
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tively difficult to observe. Few details of the spatial or 
temporal distribution or strength of their sources have 
been determined, and so they cannot adequately be de- 
scribed in global circulation models or in studies of the 
SAO. 

A likely source for high frequency gravity waves is 
convection. Alexander et al. [1995] describe obser- 
vational evidence for high frequency waves associated 
with convection, and describe the spectral properties of 
stratospheric gravity waves appearing above the deep 
convection in their simulation of a two-dimensional mid- 
latitude squall line. The stratospheric waves in their 
simulation had horizontal scales in the 10-100 km range, 
and intrinsic periods between 6 minutes and 1 hour. 
They also note that the change in direction of gravity 
wave momentum flux ahead of and behind the storm 
center may serve as a signature of waves generated by 
this source (Figure 1). This signature should be present 
in the flux carried by plane waves close to a concen- 
trated source such as a convective updraft because for 
upward conservative wave propagation, the ratio of mo- 
mentum flux to horizontal group velocity relative to the 
source is positive if the shear is not too strong. 

Pfister et al. [1993] describe some of the advantages 
in using ER-2 data for observing mesoscale waves gen- 
erated by convection. In this paper, we develop a statis- 
tical approach to the analysis of data from the Meteoro- 
logical Measurement System (MMS) [Scott et al., 1990] 
to identify features of waves with scales <100 km as an 
ensemble. Waves of this scale demand a statistical en- 
semble approach because they tend to be individually 
smaller in amplitude and in wave packet size, yet are 
ubiquitous, and so collectively may have profound ef- 
fects on the middle atmosphere. We describe a method 
for estimating momentum flux carried by 10-100 km 
scale waves along the ER-2 flight track, and compare 
the results to the Alexander et al. [1995] simulation. A 
further comparison of ozone and temperature measure- 
ments provides a consistency check on the interpreta- 
tion of the perturbations as due to gravity wave vertical 
parcel displacements. 

2. MMS Flight Data Analysis 

2.1. Momentum Flux 

Gravity wave momentum fluxes are computed from 
MMS wind vectors (u, v, w ) after application of a 10-100 
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Figure 1. Vertical flux of horizontal momentum, 
Pqu'w' as a function of distance from storm center on 
the 75 hPa pressure surface (16 km). These fluxes are 
derived from the two-dimensional squall line simulation 
in Alexander et al. [1995], and should be similar to 
those observed along an ER-2 flight path over a strong 
convective point source. The plane waves that carry 
the flux, point to their source at x=0 via the change in 
sign from left to right of storm center. The flux has been 
smoothed for comparison to the data shown in Figure 2. 

km band-pass filter to the Fourier decomposed fields. 

The filter is a smoothed version of an ideal band-pass 
[Percival and Walden, 1993] that removes mean flow 
and large scale wind variations, including most inertio- 
gravity wave perturbations, as well as small scale fluc- 
tuations that are more likely contaminated with noise 
and turbulent air motions. Filtered winds ( u ', v' , w ') are 
then perturbations at the 10-100 km scales. From these, 
vertical flux of horizontal momentum po(u'w', v'w') is 
computed and smoothed with a 10-km running triangu- 
lar smoothing function. Note that the flux is not hor- 
izontally averaged over wavelength so phase variations 
will be evident in the result. 

Time along the flight path has been converted to dis- 
tance via the aircraft speed. The Fourier analysis as- 
sumes that the wave pattern is stationary along the 
flight leg. Flying at 200 ms -1 , the plane travels 400 
km in 33 min. Violations of the stationarity assumption 
may affect the long wavelength cutoff of the bandpass 
filter but will not affect the velocity correlations because 
the aircraft speed is much larger than the wave phase 
speeds. 

Data from Flight 9 of the STEP-Tropical experiment £ 
[Russell et al, 1993] is chosen to illustrate the tech- 
nique. The initial overpass (at 18 km, 75 hPa) crosses 
a roughly linear cloud band. The aircraft data report a 
steady flight altitude along the straight 400 km long leg 
from 11°S 129. 5°E to 8°S 139°E (see Figure 2) that 
runs perpendicular to the cloud band. Following this 
high altitude leg, the inertial navigation system failed, 
so we have no MMS data for subsequent passes over this 
cloud shield. 

The momentum flux vectors along this flight path are 
shown in Figure 2. Maximum values are 0.12 kg m -1 s -2 
and northeastward in direction at 10.4°S latitude. The 
momentum flux can be seen to point in the opposite 
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Figure 2. Momentum flux vectors, p 0 (u'w', v'w 1 ) along 
the flight path shown as a function of latitude and lon- 
gitude. Each vector begins at the point of measurement 
on the flight path (the line running south-southwest 
to north-northeast). The maximum values point away 
from a region of very active convection at 10.5°S. 

direction just south of this position. The flux points 
to the northeast throughout the 60 km north of 10.5°S, 
and to the southwest over most of the 60 km south of 
this position. This change in the direction of momen- 
tum flux occurs where the coldest, highest clouds ap- 
pear below the aircraft. The apparent rotary character 
of the flux vectors is caused by smoothing effects where 
the flux direction changes sharply, and is not likely to 
be physically meaningful. 

We have no infrared cloud brightness temperature 
measurements on the 18 km leg, however, the aircraft 
turned and retraced its path on a 15.2 km leg crossing 
the same latitude 48 min later. (Note that the mean 
wind is perpendicular to the flight track and parallel 
to the cloud band which is at an approximately con- 
stant latitude during this time.) Brightness tempera- 
ture in the 10.6 pm channel [Liou et al., 1990] of the 
15.2 km leg is shown in Figure 3. On a subsequent 16.4 
km leg, the pilot had to climb to avoid flying through 



Figure 3. Brightness temperature along the flight 
track (shown as a function of latitude) from the down- 
ward looking narrow field of view radiometer 10.6 pm 
channel. The lowest temperatures show the location of 
the coldest highest cloud tops, and indicate the region 
of deepest convection. 
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an overshooting cloud turret that appeared at 10.5°S - 
10.4°S, an indication of strong vertical motions within 
the cloud below this position, and hence a likely source 
of gravity waves. We conclude the change in direction 
of momentum flux either side of this position points to 
a source below. Further to the north, fluxes generally 
decrease but are northward in direction, pointing away 
from the highest clouds, consistent with the source for 
these wave motions in the deep convection just to the 
south. Even further north along the flight path, the root 
mean square (rms) magnitude of the flux decreases by 
a factor of 4 after the aircraft passes beyond the edge 
of the high cloud band at ~ 9.3°S. This sharp transi- 
tion is further indication that the flux is carried by high 
frequency waves which propagate along more vertically 
oriented rays than the lower frequency inertio-gravity 
waves. 

2.2. Coherency of Ozone and Temperature Mea- 
surements 

Simultaneous ozone mixing ratio (Xq 3 ) measurements 
along the flight path [Proffitt et al., 1989] show a very 
high degree of correlation with temperature (T) ob- 
servations by the MMS, a property noted by Pfister 
et al. [1986] in a similar set of measurements. The 
Flight 9 band-pass filtered data show this coherency 
quite clearly. (Compare Figure 4a with solid line in 
4b.) This high degree of coherency provides an oppor- 
tunity to test the interpretation of the observed velocity 
perturbations as due to gravity waves. Any errors in the 
aircraft altitude would in turn indicate errors in the ver- 
tical velocity and vertical momentum flux because the 
winds are determined as the small residual difference of 
two large vector quantities: the ground-relative aircraft 
velocity, and the aircraft-relative wind velocity [ Scott et 
al., 1990]. The reported altitude in these measurements 




Figure 4. Perturbations along the flight track in (a) 
ozone mixing ratio and (b) temperature, shown as a 
function of latitude. In (b) the solid line is T' measured 
by the MMS, while the dashed and dotted lines show 
the results of Eqns. 1 and 2 respectively. 


is pressure altitude, Zp = H ln(po/p), where H is the 
scale height, p is pressure, and po is surface pressure. In 
the convection simulation described in the introduction, 
variations in z p are ~ ±10-20 m along the constant geo- 
metric altitude surface in Figure 1. These are the same 
magnitude as variations in z p reported along the flight 
9 leg, so we treat geometric altitude along this leg as 
constant within ~ ±10-20 m. 

The coherency between Xo 3 and T perturbations 
could be explained as either due to vertical adiabatic 
parcel displacements or aircraft vertical motions, since 
at this altitude (18 km) both Xo 3 and T have positive 
vertical gradients in their mean state profiles. However, 
the relative amplitudes in the perturbations X' 0i and T' 
differ in these two scenarios: The parcel displacement 
model includes the effects of adiabatic heating and cool- 
ing on T', while the aircraft motion model would sense 
only the vertical gradients. 

We estimate the vertical displacement A z‘ associated 
with the observed Xq 3 via: 



where dX o 3 / dz is the gradient in the background ozone 
mixing ratio. Vertical displacements derived in this 
manner for Flight 9 are ~ ±50-200 m, an order of mag- 
nitude larger than the variations in the reported pres- 
sure altitude variations. 

Two temperature perturbation models are then com- 
puted from these A z': 

(1) Assume A z' are changes in aircraft altitude. 

T- = - A *' • d i M 

(2) Assume A z' are vertical gravity wave parcel dis- 
placements. 

r GW = —Az' ■ + £) ( 2 ) 

dT/dz is the gradient in the mean state T field, and 
g/c p is the adiabatic lapse rate. The mean state temper- 
ature gradient is taken from the microwave temperature 
profiler observations [Gary, 1989] smoothed to capture 
only the large scale variations. dT/dz varies between 1 
and 5.5 °K km -1 along the flight path. Potentiahtem- 
perature gradient dd/dz is easily obtained from dT/dz. 
The ozone vertical gradient is then assumed propor- 
tional to dO /dz via the ratio of the ozone and potential 
temperature gradients determined from the aircraft as- 
cent profiles (linear fit above 16.6 km). dXo 3 /dz varies 
between 220 and 320 ppbv km -1 along the horizontal 
flight leg. 

Comparison of the models described by Eqns. 1 and 
2 above and the observed temperature perturbations 
are shown in Figure 4b as the dashed and dotted lines, 
respectively. The vertical parcel displacement model, 
T' gw (dotted line), provides a much better fit to the 
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observed temperature perturbations, and suggests grav- 
ity wave vertical parcel displacements of ~50-200 m. 
The standard error estimate describing the difference 
between Eqn. 2 and the observed T / is 0.2°K, which is 
less than the accuracy in the observations [Scott et ai, 
1990]. Other errors in this type of analysis could arise 
from errors in the vertical ozone and temperature gra- 
dients. Accuracy in 8T jdz is reported as 10% [Gary, 
1989]. Errors in the ozone gradient are expected to be 
comparable. 

It is also noteworthy that the rms values of perturba- 
tions in both ozone and temperature decrease suddenly 
by factors of 4-5 north of ~9.3°S, just as the momentum 
flux decreased at this same location. As can be seen in 
Figure 3, this location corresponds roughly to the edge 
of the high cloud band. 

3. Discussion and Conclusions 

The observed momentum fluxes in Figure 2 quali- 
tatively compare favorably with the theoretical model 
results shown in Figure 1. We see the momentum flux 
point in opposite directions on either side of a likely 
source region. The observed fluxes are, however, smaller 
by a factor of three than those in the 2-d model. This 
might be expected for any of a number of reasons: The 
2-d model might simply overpredict amplitudes in 3-d; 
the flight path may not sample the region of highest 
gravity wave flux; shear in the real lower stratosphere 
can limit wave propagation; and partial cancelation via 
superposition of waves carrying flux of opposite sigh 
can limit observed net magnitudes. The model is a 
highly idealized, 2-d, midlatitude storm. Simulations 
with more similar conditions to the observations might 
best be compared to a statistical sample of fluxes over 
convection observed with the MMS on the ER-2 in fu- 
ture work. 

Amplitudes of perturbations in the ozone mixing ra- 
tio and temperature measured along the flight path 
agree with the interpretation of these perturbations 
as due to gravity wave vertical parcel displacements. 
This may serve as an aid in interpreting the momen- 
tum fluxes determined with the method described in 
this paper. Large discrepancies in the T / - T' aw com- 
parison in other measurements could point to potential 
errors in the determination of the aircraft altitude, and 
hence errors in the vertical velocity and vertical flux of 
horizontal momentum. 

The rms perturbations in momentum flux, ozone, and 
temperature all decrease by a factor of 4-5 when the 
plane flies' north of the edge of the high cloud band, 
compared to measurements directly over the high cloud 
band. This observation may be related to the likely 
high intrinsic frequencies of the small scale waves in- 
cluded in these observations. According to the gravity 
wave dispersion relation, such waves would propagate 


at fairly steep angles from the horizon, and not be de- 
tectable far from their source at the altitudes of the 
aircraft measurements. 
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